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Abstract 

The energies of four-quark states are calculated for geometries in which the 
quarks are situated on the corners of a series of tetrahedra and also for ge- 
ometries that correspond to gradually distorting these tetrahedra into a plane. 
The interest in tetrahedra arises because they are composed of three degen- 
erate partitions of the four quarks into two two-quark colour singlets. This 
is an extension of earlier work showing that geometries with two degenerate 
partitions (e.g. squares) experience a large binding energy. It is now found 
that even larger binding energies do not result, but that for the tetrahedra the 
ground and first excited states become degenerate in energy. The calculation 
is carried out using SU(2) for static quarks in the quenched approximation 
with j3 = 2.4 on a 16 3 x 32 lattice. The results are analysed using the corre- 
lation matrix between different euclidean times and the implications of these 
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results are discussed for a model based on two-quark potentials. 
PACS numbers: 11.15.-q, 12.38.-t, 13.75.-n, 24.85.+p 
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I. INTRODUCTION 



In many-particle systems it is often convenient, or even necessary, to replace the funda- 
mental two-particle interaction by an effective interaction, which can be very different from 
the original. For example, in metals the repulsive coulomb interaction between the valence 
electrons - because of the presence of the underlying ionic lattice - is replaced by an effective 
interaction that is attractive. Also in nuclei, the free two-nucleon potential gets strongly 
modified by the presence of the other nucleons. In both of these examples, after this removal 
of the explicit photon and meson degrees of freedom, the resultant effective interaction is still 
mainly in the form of a two-body interaction - with only a minor three-body term arising 
in the nuclear case. However, for a system of quarks interacting via gluon-exchange - even 
though this basic interaction is that of QCD - little is known in multiquark systems about 
the corresponding effective interquark interaction after the explicit gluon degrees of freedom 
have been removed. This is an important question, if realistic calculations are to be made 
for interacting quark clusters e.g. as in meson-meson, meson-nucleon or, eventually, nucleus- 
nucleus scattering. At present, the only way to carry out these fundamental calculations is 
by using Monte Carlo lattice techniques. Unfortunately, mainly due to the limitations of 
present day computers, these calculations are restricted to clusters of, at most, two or three 
quarks - see for example Ref. JTJ. Therefore, a bridge is needed between lattice calculations 
involving few quarks and conventional many-body techniques that can accomodate larger 
numbers of quarks. 

In an attempt to throw some light on the relationship between QCD and the effective 
interquark interaction, a series of model calculations have recently been made @- ||. In 
these, the energies of four quarks in various configurations (e.g. on the corners of a rect- 
angle or tetrahedron) have been calculated in quenched static SU(2) on a 16 3 x 32 lattice. 
The reason for studying four quarks is because of the possibility for partitioning these into 
different colour-singlet groups each containing two quarks - a situation not possible with 
2 and 3 quark systems. This can then be considered as a step towards the scattering of 
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quark clusters - in this case meson-meson scattering. Hopefully, many of the features of 
an effective interquark interaction in SU(2) will be preserved in the more realistic case of 
SU(3) - in the same way as, for example, the ratios of glue-ball masses/string energy are 



numerically very similar in SU(2) and SU(3) ||- ||11|| . 

In Refs. ||- ||, the main feature in the results is that the strongest interaction between 
two separate two-quark clusters occurs when two of the three possible partitions of the 
clusters are degenerate, or almost degenerate, in energy. This was first observed in Ref. 
H], where the four-quark binding energy E dropped by an order of magnitude, when the 
four-quark geometry changed from that of a square [r x r] to a neighbouring rectangle 
[r x (r ± a)} - where r is the length of a side of the square in lattice units a. For example, 
E(3a x 3a) = -0.054(l)/a = -90 MeV compared with E(3a x 2a or 4a) « -0.006(l)/a. 
Subsequent works verified this observation with other geometries, where the two partitions 
were not exactly degenerate - but more so than for the above rectangles. For example, in Ref. 
this was achieved by tilting the rectangles out of the planes defining the underlying lattice, 
and also with other non-planar geometries. These latter cases showed that it was indeed the 
energy degeneracy dictating the size of the interaction and not a geometrical degeneracy. 
In other words, with clusters (13) (24) and interquark potentials V(ij), it is the degeneracy 
V(13) + V(2A) « V(U) + V(23) that matters and not |r x -r 3 | + |ri -r 3 | « |r x — r 4 | + [r 2 — r 3 |. 
For large interquark distances V(ij) ~ & s |r.; — Tj\ + c, where b s is the string energy and c 
an additive constant. In that case, the energy- and geometrical- degeneracy are the same. 
However, this is only true for > 0.5 fm and in the interesting region dominated by 
explicit quark interactions it is the energy degeneracy criterion that must be used. This is 



contrary to the assumption made with the so-called flip-flop model of Ref. |12j that takes 
the geometrical degeneracy criterion for all interquark distances. 

The reason for studying many different 4-quark configurations (i.e. rectangles, linear, 
non-planar, etc.) is to give a representative selection of geometrical possibilities that arise 
in practice when, for example, two mesons scatter. Any quark model for such a scattering, 
when reduced to the same conditions as the lattice calculation [i.e. static, quenched, SU(2)] 



should give for the above selection of configurations the same values for the 4-quark energies 
E. As said earlier, the representative selection of 4-quark configurations concentrated on 
those cases where two of the three possible partitions were degenerate in energy. However, 
a potentially interesting case is that of the tetrahedron, where all three partitions are now 
degenerate. This is the main subject of this paper and is discussed in section 2. For 
completeness, the series of configurations shown in figure 1 are treated with 1 < d < 5 and 
1 < r < 6 - the tetrahedron being r = d. The r = 1 case is then closely related to the 
earlier work on squares with d = 1, . . . , 5 and serves as a check on those calculations, since 
they were done with a set of basis states restricted to two in number and not three as in the 
present work. 

Since the completion of Refs. ||- J?|, it has been suggested |L3| that perhaps more care 



is needed in extracting energies and their statistical errors from a series of Monte Carlo 
results. Therefore, in section 3, the data is analysed taking into account data correlations 
between different euclidean times. 

II. THE LATTICE MONTE CARLO CALCULATION 



Since most of the details of the lattice Monte Carlo calculation have been given in the 
earlier papers 0- |J, this section will only very briefly outline those details concentrating 
mainly on aspects specific to the geometry shown in Fig. [I]. 

Throughout this study the four quarks are treated in the static quenched approximation 
on a 16 3 x 32 lattice in SU(2) with f3 = 2.4 - corresponding to a lattice spacing of a ~ 0.12 fm. 
Both the four-quark energies V^(4g) and the corresponding two-quark energies Vi(2q) are 
extracted on this lattice. However, the main quantities of interest are the ground state 
binding energy E i=1 and the excited state energies E i=2 ^ defined by the differences 

Ei = Vi(4q) - 2Vt(2q) (1) 

As described in Refs. ||- f5[, the V{2q) potentials are calculated from a 3*3 variational 
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basis, where - after fixing the gauge so that all temporal links are set to unity - the three 
states are constructed from basic lattice links that have been "fuzzed" to different levels -12, 
16 and 20 in the spatial directions. This fuzzing essentially models a glue flux-tube between 
the two quarks in question. For the cases where these two quarks are not along a given 
spatial axis [e.g. V(x, y)], the appropriate path of links is then constructed as the average of 
the two most simple paths connecting x and y - each consisting of one straight section along 
the x and y axes. For V(4q), using only the maximum fuzzing level of 20, the basis states 
are constructed from the different partitions of four quarks into a product of two two-quark 
colour-singlet clusters. In some cases it seems natural to restrict the variational basis to just 
two partitions e.g. with rectangles the two involving the sides of the rectangle and not the 
one constructed from the two diagonal paths. However, in other cases - especially for the 
tetrahedron of interest here - it is seems to be necessary to keep all three partitions. 

Having constructed the above paths (e.g. PiiTi) = t/i(Ti)|Vac) and -Py(T 2 ) = 
Uj(T 2 ) | Vac)) from fuzzy links, Wilson loops - the quantities from which the final ener- 
gies are extracted - are then simply given, with the present gauge, as the overlap of paths 
separated by euclidean time T = T 2 — Ti i.e. = (Pj(T\ + T)\Pi(Ti)) . 

The energies Ei in Eq. (|TJ) are expected to be exact in the sense that they are not 
the result of a truncated weak or strong coupling expansion. However, they do contain 
uncertainties such as: 

a) Statistical errors. To a large extent these can be estimated with reasonable accuracy 
- the subject of section 3. 

b) Systematic errors. By their very nature they are much harder to estimate. Through- 
out, the aim is to attempt at minimizing their effect. In general, the most obvious source 
of such errors is the use of a lattice that is too small (i.e. too few points) or has a grid that 
is too coarse (i.e. (3 too small) for the results to represent those expected in the continuum 
limit of the lattice. This was checked in Ref. || by using a 24 3 x 32 lattice at /3 — 2.5, which 
corresponds to a lattice spacing a ~ 0.082 fm. In the present type of calculation, since - 
as seen from Eq. (|l|) - the signal E{ is a delicate cancellation, it is necessary to ensure that 
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the appropriate Vi(4q) and Vi(2q) are extracted simultaneously from the same lattice 
configuration. In this way, systematic errors possibly present in the separate energies Vi(4q) 
and Vi(2q) are expected to cancel to some extent. 

In order to study different aspects of degenerate geometries, three separate sets of lattice 
calculations are performed in this paper. 

1 . The first concentrates on geometries near to the tetrahedra - namely - those cases 
with r — d, c?± 1 (except for d = 1, where only r = 1, 2 are treated). The data is 
collected into 69 blocks each containing 32 measurements i.e. 2208 measurements 
in all. 

2. The second set of geometries is for r = 1 and d = 2 to 5 i.e. situations close to 
that of the squares studied earlier in Refs. 0- [[/] . Here the data is collected into 
47 blocks each containing 64 measurements i.e. 3008 measurements in all. 

3. In order to see the connection between the present series of geometries - especially 
that in the previous item - and the square/rectangle geometry of Refs. [[J- ]7j|, a 
third set of runs is carried out. These are a repeat of the earlier square/rectangle 
geometry but with the complete 3*3 basis i.e. with the inclusion of the basis state 
constructed from the diagonal paths. Only a limited number of 64 measurements 
are made - 4 blocks each containing 16 measurements. However, this is sufficient 
to see how the energies vary in going from the three 2*2 bases (A+B, A+C and 
B+C) to the complete A+B+C basis. 



III. ANALYSIS OF THE DATA 



As emphasized, for example, in Refs. ||13[| - fH}| , care should be taken when analyzing data 



generated in Monte Carlo lattice calculations, in case the data between different euclidean 
times is too correlated. This can be taken into account in the following manner. As discussed 
in the previous section, the actual quantities measured are the Wilson loops between 
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different states i,j at different time intervals T. Here the i,j refer to the different degrees 
of fuzzing for the two-quark potentials and, for the four-quark case, to the different quark 
partitions. However, the quantities of interest are the energies (V) of these systems. In Refs. 
||- |7|] these were extracted by solving the eigenvalue equation 

W?aJ = X^W^af, (2) 

where A- T ^ — > exp(— Vi) as T — > oo. This matrix formulation takes into account mixing 
between the different degrees of fuzzing or four-quark partitions, but says nothing about 
possible correlations between the different times T. As suggested in Ref. |I3J it is convenient 
to treat the configuration and time correlations in the same manner. For clarity, this is 
illustrated by the specific example of n partitions and T time differences in the four-quark 
case. This involves the (n 2 T) matrix elements W^, W± 2 , ■ ■ ■ , W% n , where the expected equal- 
ity between, for example, W± 2 and is not enforced by symmetrizing the data. Let the 
index p = 1, . . . , n 2 denote these n 2 matrices. As discussed in the previous section, the data 
(W?) in the form of these (n 2 T) matrix elements is collected into N separate blocks, each 
of which is already the average of m measurements i.e. in all there are mN measurements 
of each Wilson loop. The average of these measurements is denoted by 

iV i=i 

The problem is now reduced to fitting these averages by a function of the general form 

k M 

if = exp(-V fc T) (4) 

fc=i 



in order to extract the energies Vk- As shown in Ref. |13[ this can be achieved by minimizing 
the expression 

X 2 = E {F^-W^M(T,T',p,p')(F^-W}'), (5) 



T,T>,p,p> 



where 



M(T,T',P,P') = NC-\T,r,p,p f ) (6) 



and C is the covariance matrix 

N 



C(T,T>,p,p>) ~EKW - «(0 - O • (7) 



=1 

In what follows, it is convenient to introduce the correlation matrix 



Crr.r. „.,/)■. c(T,r,p,j/) 

IC{T,T,p,p)C{T',T>,p,p) 



The above procedure is stable provided A" is large enough. However, in the present type 
of lattice calculation this is not guaranteed and can then lead to C(T,T',p,p') having very 



small eigenvalues that produce large fluctuations in M(T,T',p,p'). As argued in Ref. |T3 



such eigenvalues correspond to eigenvectors that alternate in sign as a function of T and so 
are not very relevant to smooth fit functions of the type proposed in Eq. (ED). It is, therefore, 



reasonable to remove these disturbing eigenvalues and here the suggestion made in Ref. [13 
is adopted. This simply replaces the smallest eigenvalues of C(T,T',p,p r ) by their average. 
The question then arises concerning how many (n t ) of the original eigenvalues = n 2 T) 
should be retained. Initially a number tim of these Ut eigenvalues are fixed - with um being 



chosen as ~ v A" - a compromise choice found after some experimentation in Ref. [|T3j . The 
average (X av ) of the remaining n T — n M eigenvalues is then made and those eigenvalues less 
than (A a „) are replaced by (X av ) itself. Since only a fraction of the n T — n M eigenvalues are 
replaced in this way, there still remains n t > Um of the original eigenvalues. In practice 
the final value of n t is not strongly dependent on the original choice of um- In the present 
problem A" = 69, n = 2 or 3 and T can range from 1 to 5, so that for the largest C matrix 
only about 8 of the total number of 45 original eigenvalues should be retained. However, 
frequently the T = 1 data are dropped from the analysis, since they have very small error 
bars and presumably contain the effect of several higher energy components that become 
negligible already at T = 2 and so cannot be seen for T > 2. Therefore, keeping the T = 1 
data often makes it difficult to obtain an accurate fit with only a few (1-3) terms V^. 

Having arrived at a suitable model for C, it remains to fix the precise form of the function 
Fp in Eq. (RJ), which is to fit the n 2 T matrix elements by means of Eq. (|5|). It is seen that 



the Fp are described by the ku potentials and the corresponding r^ku amplitudes a p (k) 
i.e. in all (n 2 + 1)/cm parameters. For the four-quark case the most reliable fits are found 
with three partitions (n = 3) and four time steps (T = 2, . . . , 5), which means 36 pieces of 
data are to be fitted with 10/cm parameters. Therefore, in principle only three potentials 
(^1,2,3) can at most be extracted. However, fitting 36 numbers (some having significant error 
bars) with 30 parameters (&m = 3) would not succeed. In fact, even the extraction of two 
potentials using the 20 parameters (&m = 2) would probably not give values of Vi,2 with a 
meaningful accuracy. It is, therefore, necessary to impose some theoretical constraints to 
restrict the number (n 2 kM) of parameters a p (k). Since the data Wf^ involve the overlap of 
two lattice configurations Pj(0) = C/j(0)|Vac) and Pj(T) = Uj(T)\V&c) separated by time T, 
it can be written as the vacuum expectation 

W?j = (Vac|C/ j (T)C/ i (0)|Vac) = ^(Vac^T) \k)(k\Ui(0)\ Vac) exp(-V k T) (9) 

fc 

It is, therefore, reasonable to assume that the amplitudes are separable i.e. a p (k) = 
a ij(k) = di(k)aj(k). This has two good features: 

1) The number of parameters is now reduced to (n + l)k M . Therefore, fitting the 36 
pieces of data from three partitions (n = 3) and four time steps now only requires Aku 
parameters. This makes the quite reliable extraction of two or possibly three potentials 
feasible. 

2) This parametrisation imposes the symmetry Fj- = Fj t as is expected physically. As 
said earlier this symmetry was not forced on the original data by some procedure such as 
taking the average 0.5(W? + WJ-), since the differences between and Wj- can contribute 
to the error analysis of the extracted energies. 

It should be added that, in the event of two energy states being degenerate (e.g. V kl = 
Vfc 2 ), then ai(ki)aj(ki)+ai(k2)aj(k2) should be replaced by a single amplitude a^. This is the 
situation encountered in the tetrahedron geometry. However, there additional symmetries 
suggest the following better model. 

In some symmetrical cases, such as four quarks on the corners of a square or tetrahedron, 
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the number of parameters can be further reduced by choosing forms of FT that guarantee 
various symmetries. For the square, when only the two partitions involving the sides are 
considered, the matrix of Wilson loops has the form 



( ws wl ^ 



V 



(10) 



wl wl 2/ 

where not only is the general symmetry W\ 2 = W 2 i expected but also for a square W\\ = W 22 . 
In this case Eq. (0) is easily solved to give for the lowest energy 

wl + W? 2 



Ai 



w 



T-l 
12 



(11) 



and for the energy of the first excited state 



A 9 



wir 1 - wl- 1 



(12) 



Therefore, in fitting the data it is reasonable to expect that only two potentials can be 
extracted i.e. ku = 2, and so the fitting functions must have the form 



if 



Fj = a(l) exp(-KT) + a(2) exp(-V 2 T) 



FT = FT 



a(l) exp(-KT) - a(2) exp(-V r 2 T). 



(13) 



This reduces the number of parameters to four [Vi,2, 2)] in order to fit the 16 pieces of 
data covering four T steps. Here it should be remembered that the suffix p on FT is fixed 
as p = j + n{i — 1). 

Similarly, for the tetrahedron the matrix of Wilson loops has the form 



T 



w 



where the general symmetries WT X 



WT 2 Wl 



\ 



13 



Wl Wg Wl 

wl wl wl 



wl = wl, wl = wl, wl 



wl, wl 



32 



(14) 



wl 



are expected and, in addition, there are the equalities Wl = Wl and Wl = — Wl. The 



2.3 



13- 
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" minus" sign appearing in the last equation is a reminder that the quarks are in fact fermions 
even though quarks and antiquarks transform in the same way under SU(2). This point is 
discussed in more detail in the appendix of Ref. J7|. Therefore, in all, there are only two 
independent Wilson loops and W^ 2 - Again Eq. (g) is easily solved to give for the lowest 
energy (occurring twice) 

Al,2 

and for the excited state 

A 3 = 

In fitting the data it is reasonable to expect that only two terms can be extracted, i.e. 
ku = 2, and so the fitting functions must have the form 

Ff = F 5 T = Fg = a(l) exp(-y 1 T) + a(2) exp{-V 2 T) 

F 2 T = F 3 T = Fj = Fj = —Fq = —Fg = 0.5a(l) exp(-^T) - a(2) ex P (-V 2 T). (17) 

This again reduces the number of parameters to four [Vi i2 , a(l, 2)] in order to fit the 12 pieces 
of data (n = 3) covering four T steps. In this case, due to the high symmetry there are no 
longer 36 pieces of independent data. 

In the tetrahedral case it is also seen that, when one of the partitions is removed (i.e. n 
drops from 3 to 2), the Wilson loop matrix reduces to the same form as that in Eq. (|T0|). As 
seen by comparing Eqs. (|ll]) and ([15]) the lowest eigenvalue is now the same in the n = 2 
and 3 cases. However, the first excited state is quite different in the two cases - unlike the 
result for most other geometries. For example, as will be seen later in Table [IV] of Section 
P"V| , for squares and for those rectangles so far discussed, the n = 2 and 3 cases give very 
similar results for E 2 as well as E\. However, there the third partition has a higher energy 
than the other two partitions, whereas for the tetrahedron all three partitions have the same 
energy. When the tetrahedron is distorted into a neighbouring lattice configuration, it will 

12 



Wl + W T 



12 



wTi~ 



1 + wj 2 



wl - ml 



w 



T-l 



2W, 



T-l ■ 
12 



(15) 



(16) 



be seen in Table [H] of Section |V| that both E\ and E2 can again be given quite accurately 



by n=2. It should be added that for very elongated rectangles the two partitions with the 
highest energy become more and more degenerate. In that case E\ is essentially zero and 
E2 will depend more and more on the presence of the third partition. 

In an idealised situation, all of the available n 2 pieces of data from T=l to 5 should - 
using the above procedure - be fitted with a containing ku ~ n potentials. This would 
then be the natural extension of Eq. (§) for incorporating into its n eigenvalues the effects 
of correlations between different values of T. As it now stands that equation only includes 
the mixing between the different fuzzing levels or partitions. However, in practice, fitting all 
the T values is not possible and decisions have to be made concerning both the data (W?) 
to actually be fitted and, in the fitting procedure, the form of F^ and also the model for 
the correlation matrix in Eq. (j8|). 

a) In the present problem, the data goes from T—l to 5 and the number (n) of fuzzings 
(2-quark) or partitions (4-quark) can range from 1 to 3. For both the 2- and 4- quark 
systems usually T=2 to 5 turns out to be the most suitable range. When the T range is 
reduced to 3-5, better values of \ 2 emerge from Eq. (|5|) but at the expense of larger errors 
on the extracted potentials V$. On the other hand, when T = 1 data is included - because 
of the small errors on this new data - it usually results in either a "no-solution" situation or 
a x 2 that is too large to be meaningful. In addition to the restriction on T, a decision must 
be made on the number of fuzzings or partitions to be included. In the 2-quark system only 
two fuzzings were necessary at a given time. This is simply a reflection that the different 
degrees of fuzzing are effectively very similar and so little is gained by increasing the number 
of different fuzzings from 2 to 3. Here the two largest fuzzing levels (16 and 20) are used i.e. 
n is always taken to be two in the 2-quark case. In the 4-quark case it is possible to use all 
of the data from the different partitions i.e. n is always three in the tetrahedron case. 

b) In FT , once the separable form or the more symmetrical one in Eqs. (|T3|) or (|T7|) have 



v 

been chosen, the only decision to be made is the value of &m _ the number of terms (Vj) to 
be determined. This is varied between 1 and 3. In the correlation matrix C(T,T' ,p,p') of 
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Eq. (^|) the number of eigenvalues to be initially retained (um) can also be varied. Ref. fl3 



suggests using um ~ \N, and this was found to be a suitable value in the four-quark case 
as can be seen from Table | (there iV = 69, so that the -y/iV-result is shown as case 11). On 
the other hand, in the two-quark system it is necessary to use um = 2 or 3 to get a stable 
fit. Table | shows the effect of varying % in a four-quark system (cases 9-13), as well as a 
fit with no correlations (i.e. (7=1, case 8). A few general features can now be seen: 

1. Comparing cases 8 and 11, the use of the correlation matrix has approximately a 
one-sigma- level effect on the averages and it slightly decreases the errors. 

2. The exact value of hm is not important as long as approximately all the stable eigen- 
values are retained (i.e. nt is in the right range). To estimate this range for n t a 



jackknife procedure on the eigenvalues was carried out [15 



The above choices in the data and fitting procedure are strongly inter-related and so a 
strategy is necessary to single out the optimal set of values for n, the range of T, ku and 
Um- An example of this is given in Table |. Basically this strategy amounts to first fitting 
only a portion of the data e.g. T =(3 or 4) - 5, with only a few potentials e.g. ku =1 or 
2 - see cases 1-2 in Table |. It is possible to see a suitable range for n t , and therefore for 
Um as well, already from these fits. Then ku is increased to find a maximum number of 
potentials that can be extracted by using this range of T (see cases 2,6,7). The data base is 
then progressively enlarged to include more timesteps, while trying to keep the value of x 2 
per degree of freedom fixed or decreasing by adding more potentials if possible (cases 2-7). 
Naturally, adding a timestep constrains the potentials more than before, so that the errors of 
the potentials decrease accordingly. Initially this works well, but eventually the \ 2 become 
larger and larger, so that usually the results become meaningless before the ultimate stage 
of including all of the data (i.e. n = 3 and T =1 - 5). However, it should be added that the 
T = 1 data fitted with ku = 3 has to be treated with caution, if the resulting potentials 
are significantly different from those extracted with T = 2 — 5 and A;jv/=2, since this may 
indicate that even higher energy states are polluting the T = 1 data. 
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IV. RESULTS 



The results for the tetrahedron-like geometry in Fig. [I] are shown in Fig. |2] and Tables [Tl 
and [TTT]. Several points arise from Table fj: 

1. From earlier work in Refs. [[3]- ||, for the corresponding squares [i.e. (d, 0) in 
this notation], where two of the basis states are degenerate, the binding energy 
of the lowest state ranges from -0.07 to -0.05 as d goes from 1 to 5 [see Table 



rV| . However, now - even though at least two of the basis states are always 
degenerate (A and B) - the ground state binding energy is always less than that 
of the corresponding square [E 1 (d,0)}. For a fixed d, \Ei(d,r)\ decreases as r 
increases. Nothing interesting happens to E\ at r = d, at which point all the 
basis states are degenerate in energy. 

2. For fixed d, as r increases from to d, the energy of the first excited state E 2 
decreases until E2(d,d) = Ei(d,d). For r > d, E2(d,r) increases again. This 
degeneracy of E\^ for the tetrahedron is a new feature compared with earlier 
geometries. As will emerge in the next section, this is a severe constraint on any 
model wishing to describe this data. 

3. As seen from columns A + B + C (I and II), in all cases the energies extracted 
using either (J) Eqs. (§])— © or (//) Eq. (||) are within the error bars of each 
other. If a choice has to be made, then the results (/) are prefered, since they 
involve a more detailed analysis of the data. 

4. In columns A+B and B+C only the corresponding two basis states are used. The 
symbol 5* means that the extracted E\ are exactly the same as the A-\-B + C(II) 
results. In those cases it is, therefore, sufficient to only use a 2*2 basis. However, 
the choice of which 2*2 basis depends on the particular geometry, since one of 
these two basis states must have the lowest unperturbed energy. Since A and B 
are degenerate in energy, for a given d this amounts to using A + B for r < d. 
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Whereas, for r > d it is necessary to use B + C, since C now has the lowest 
unperturbed energy. 

5. Except for the tetrahedra, the values of E 2 are essentially the same in the 2*2 
and 3*3 bases. 

6. The values of E 3 (d,r) are always much higher than E 2 (d,r). However, as dis- 
cussed in the next section, this second excited state is dominated by excitations 
of the gluon field and so is outside the scope of the models introduced in that 
section. 

7. The last column shows the two 2-quark ground state potentials Vi(d,d) and 
V\_(d,r) that enter with this geometry. In the next section, these are needed 
for evaluating a model for the 4-quark binding energies. 

Several of the above points are also illustrated in Fig. ^| for the case of d/a = 3 and 
r/a — 0, 1, 2, 3, 4. In addition, there are shown the theoretical predictions for the / = 1 limit 
of the model to be discussed in the next section. 

Table |TTT| is very similar to Table |TJ except that r/a is now restricted to unity i.e. geome- 
tries nearest to the (d x d) squares. The following points can be made: 

1. Again the energies extracted by Eqs. @— © (/) and Eq. fl2|) (II) are within 
the error bars of each other - with (/) being preferable, since it involves a more 
detailed analysis of the data. 

2. For comparison in Ei(d x d) are given the energies of the nearby squares from 
Ref. |J with only a 2*2 basis but analysed with method (I). In general, these 
all have a somewhat larger binding energy, indicating that the binding energy is 
maximized for the square geometry. 

In Table |TV| results are given for various square/rectangle geometries. The main purpose 
of this set of runs is to see any change in the results for E 1}2 extracted from the 2*2 bases, 
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when this is extended to a 3*3 basis by including a third state C. Also an estimate can now 
be made for the position of the second excited state E 3 . 
Several comments should be made about these results 

1. For the lowest energy Ei, introducing the third basis state has little effect on the 
results given by the two basis state calculation provided the partition with the 
lowest energy is one of those two states in that basis. In fact, for squares (i.e. 
r = d) the results are identical for the A + B + C and A + B bases. However, in 
a general case when using a 2*2 basis, for d > r it is necessary to include state 
A and for d < r state B. Remember that in this notation it is state C that is 
constructed from the diagonal two-quark colour singlets and so always has the 
highest unperturbed energy. 

2. The energy Ei of the first excited state is less dependent on the choice of basis 
- with the 3*3 and 2*2 possibilities giving in most cases essentially the same 
results. 

3. From column A + B + C it is seen that E 3 - the energy of the second excited 
state - is always considerably higher in energy than E 2 . 

4. The Ei(P) column is the result of many more measurements than the following 
columns. Therefore, in view of the facts in the previous items, these numbers 
should be considered as the most accurate for rectangles. They are basically the 
same as those in refs. ||, and 0, except that the data has been subjected to 
the improved analysis (J) using Eqs. ©-(H). 

V. AN INTERPRETATION OF THE RESULTS 

One of the main reasons for embarking on this work is the attempt to find a model which 
gives a simple understanding of the four-quark binding energies in terms of the corresponding 
two-quark potentials. At first sight, it may seem that the best such model should be able 
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to explain all of the energies extracted from the lattice calculation of the previous section. 
However, this would be not only too ambitious but also it would be outside of the goal of 
the desired model. When the lattice data is analysed by directly diagonalizing Eq. (fj), then 
the number of energies extracted is equal to the number of basis states i.e. three for the 
tetrahedral geometry. On the other hand, when the expression in Eq. (|5]) is minimized, the 
number of extracted energies depends very much on the quality of the lattice data. For the 
tetrahedron, as seen in Tables p lTJ, it is also possible to extract the three lowest eigenvalues 
and, in principle, even more could be obtained. This point is discussed in more detail in ref. 
fl6|| , where also the results from other four-quark geometries are extracted using Eqs. (|5])- 
(|). However, for the model to be discussed below only the lowest two lattice eigenvalues 
are of interest, since the third is dominated by excited gluon components. Such a feature is 
outside any model that incorporates an interaction that is based on the two-quark potential 
in its ground state. It can be seen that the third state is basically a gluonic field excitation 
as follows: 

1) As shown in the above tables and also in the appendix of Ref. J7|], except for the 
tetrahedra, the addition of a third basis state into the lattice calculation has only a minor 
effect on the results using only two basis states - indicating that the structure of the third 
eigenstate is very different from the lower two states. 

2) In the linear case, on a lattice where all three partitions are constructed from links 
along a single spatial axis, i.e. without the need for introducing combinations of links along 
different axes, these partitions are linearly dependent leading to a singular matrix, if a three 
basis calculation is attempted in Eq. (0). 

The previous paragraph has shown that only the lowest two lattice eigenvalues are of 
interest when constructing a model based on the basic two-quark potential in its ground 
state. To this end, in refs. 0- it was proposed that the lattice energies E^ should be 
fitted by a model defined in terms of only two basis states A and B. In this model the Ei 
are given as the eigenvalues of the equation 
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[V(/) - Ai(/)N(/)] % = 0, 



(18) 



with 



N(/) 



1 //2 
//2 1 



and V(/) 



^ «13 + ^24 fV AB 



V 



(19) 



where A« = Ei + v± 3 + f 24 - the configuration A = {qiqz){q 2 qi) being the lowest in energy of 
the three possible partitions into two two-quark singlets. The off-diagonal matrix element 



(A\V\B) = V AB = V BA = - (u 13 + v 24 + v 14 + v 23 - v 12 - v 3A ) 
is of the form expected from a two-quark isovector potential 



(20) 



Vij — Tj • TjVij. 



(21) 



The additional factor of / appearing in the off-diagonal matrix elements is to be interpreted 
as a gluon field overlap factor. In the weak coupling limit / is unity. However, in general this 
limit is found to result in too much binding compared with the lattice results. Therefore, 
/ is treated as a phenomenological factor, which is adjusted to fit the lattice energies. The 
problem is then reduced to understanding the resulting values of /, which have essentially 
simply replaced the lattice energies. That this is a reasonable model is supported by several 
points: 

1. This factor / is approximately unity when all four quarks are close together. This 
is indeed where the f = 1 weak coupling limit should be best. 

2. For a given geometry, a single value of / qualitatively explains quite well both 
Ei and E 2 . 

3. When / is parametrized as 



/ = exp[-kb 3 S\, 



(22) 



where S is the minimal area of a surface bounded by the straight lines connecting 
the four quarks and b s a 2 is the string energy density with the value ~ 0.070, 
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then it is found that - for squares and rectangles - the parameter k is reason- 
ably constant at w 0.50 ± 0.05. Of course, the ultimate goal is to find such 
a parametrization in which the corresponding parameter(s) k would be strictly 
constant for all geometries (i.e. squares, tetrahedra,. . .). If this selection of ge- 
ometries is sufficiently representative, then the reasonable assumption would be 
that this same parametrization should work for other geometries not calculated 
on the lattice. 



Prior to this work on tetrahedrons the geometries considered had, at most, two of the three 
possible partitions being degenerate in energy (e.g. for squares). In these cases, as seen from 
Table [IV], it is found that the lattice energies E\ and Ei are essentially the same for the 
three-basis-state calculation (A+B + C) and those two-basis-state calculations (A+B, A+C 
and effectively B + C) which involve the basis state with lowest unperturbed energy. This is 
one of the reasons why the 2*2 version of the /—model in Eq. (|l^) was quite successful for 
a qualitative understanding of these cases. However, for tetrahedra and the neighbouring 
geometries calculated in Section 2, it now seems plausible to extend the /-model to the 
corresponding 3*3 version in which 



N(/) 



1 f/2 f/2 
f/2 1 -f/2 
[ f/2 - f/2 1 j 



^ V 13 + V 2 4 fV AB f'V AC ^ 



and V(/) 



V 



fV B A v u + V 2 3 -f"V B c 
f'VcA -fVcB V 12 + f 3 4 



(23) 



where the negative sign in the BC matrix elements is of the same origin as the one in Eq. 



(|14]). This extension has both good and bad features. On the positive side, all three basis 
states are now treated on an equal footing. This is convenient when considering some general 
four-quark geometry, since it is then not necessary to choose some favoured 2*2 basis, which 
could well change as the geometry develops from one form to another. On the negative side, 
in the weak coupling limit (i.e. /, /', /" — > 1) the 3*3 matrix of Eq. (|T8| ) becomes singular 
- in the sense that adding columns B and C results in column A. However, in this limit, 
each of 2*2 matrices corresponding to the three possible partitions A+B, A+C and B+C 
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gives the same results. Away from weak coupling the 3*3 matrix is no longer singular, but 
now the three possible 2*2 partitions do not necessarily give the same results. Below an 
attempt is made to minimize the differences between these three partitions, since in the 
corresponding lattice calculation the differences in most cases are indeed small. 

The strategy of trying to mimic the lattice calculation by means of the /-model in Eqs. 
([[ID and (|19D is carried out at two different levels. 

1. The lattice calculation is made in the static quenched approximation with the 
SU(2) gauge group. Therefore, in the /-model there must not be any kinetic en- 
ergy term (i.e. static quarks) and also quark-antiquark pair creation (for example 
in the form of meson exchange between quark clusters) must not be included (i.e. 
the quenched approximation). Furthermore, the two-quark potential in Eq. fl21~[) 
is expressed in terms of r-matrices - the generators of SU(2). 

2. The lattice calculation, as said earlier, gives essentially the same results for any 
of the three partitions A+B, A+C and B+C as the complete A+B+C basis. As 
will be seen below, this point is more difficult to mimic. Since this feature is 
automatically encoded in the matrix elements of Eq. ([14]), one possibility for the 
/-model is to attempt to mimic in more detail the form of these matrix elements. 



For the tetrahedron geometry, the /-model as written in Eq. (23) takes on a particularly 



simple form since each of the are equal and also / = /' = /". In this case the eigenvalues 



are 



E 1 = E 2 = E 3 = 0. (24) 

This has the positive feature that E\ and E 2 are degenerate as in the lattice results of Table 
[1]]. However, they are degenerate at zero binding energy. This latter feature is unavoidable 
for the model in its present form, since there is only one energy scale present. This is the 
two-quark potential between each of the quarks, which is, of course, independent of the 
pair of quarks chosen, since all interquark distances are the same for the tetrahedron. It 
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is, therefore, necessary to introduce a second energy scale into the model. However, any 
improvements in the model have very limited choices, since there are only two different 
matrix elements involved - the diagonal ones all equal to — E and the off-diagonals ones 
all equal to ±0.5 fE. Therefore, the most general modifications are to change the diagonal 
matrix elements to d\ — E and the off-diagonal ones to ±0.5/(cf 2 — E). This results in the 
eigenvalues 

di + 0.5fd 2 rfi - fd 2 , , 

E 1 = E 2 = 1 + Q5f and E 3 = -j-^. (25) 

At first sight it may appear that there is sufficient information to now extract the new 
parameters dip, since / can be estimated using the parameters (assumed to be universal) 
from other geometries - thus leaving two equations for Ei 2 and E 3 and the two unknowns 
dip. However, as said before, this is too much to demand from the /-model, since in 
the lattice calculation the third basis state in the complete A+B+C basis generally plays a 
minor role in determining the values of Eip and, therefore, the third eigenvalue is presumably 
dominated by an excitation of the gluon field. Even so, it is of interest to see that a similar 
feature now arises with E 3 in Eq. ([25]) , since this third state is removed in the weak coupling 
limit i.e. E 3 — > oo as / — > 1. However, the /-model - if it is to be successful - should 
only be expressed in terms of the lowest energy gluon configurations, since the gluon field is 
not explicitly in its formulation, but only appears implicitly in the form of the two-quark 
potentials and the /-factors. In view of this, no quantitative attempt should be made to 
identify the second excited state emerging from the lattice calculation as E 3 in the /-model. 
In Ref. JT/J it was shown that the two-state model of Eqs. (|T9|) with the overlap factor 



/ = 1 agreed with perturbation theory upto fourth order in the quark-gluon coupling [i.e. 
to 0(a 2 )} and gives Ei 2 =0 for tetrahedra. Therefore, the non-zero lattice results for small 
tetrahedra must be of 0(a 3 ) at least. Another aspect of this special situation for tetrahedra 
is also seen - when extracting or interpreting the value of E 2 - by the need for the third basis 
state both in the lattice calculation (see Table [Tip and in the /-model, since in comparison 
with Eq. ( P5|) the two basis state version gives 
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= 4+0.5/4, md = j, - 0.5/4, 

1 + 0.5/ 1-0.5/ V ; 

i.e. both the two- and three- basis state models have the same ground state, but the latter 
does not show the E\ = E 2 degeneracy. 

In Fig. |2] the predictions are shown for the / = 1 limit of the above model in the case 
of d = 3. If, in the notation of Fig. [I], the appropriate two potentials in state A(or B) and 
state C are defined as v a = Vi(3,r) and v c = Vi(3,3), respectively, then 

Vl{M) = 8v * ~ 2v c and V 2 {M) = 2v c . (27) 

Here the V^(4g) are defined in Eq. ([!]). To extract Ei for r < d the two-quark potential 
Vi(2q) is taken to be v a , whereas for r > d the appropriate potential is Vi(2q) = v c . 

The expressions in Eq. (|25|) are not particularly useful unless there is a model for the 
parameters dip- However, since it's not the purpose of this paper to make a comprehensive 
study of models covering all the 4-quark geometries considered in earlier works [|J- ||, only 
a few general remarks will be made here for the tetrahedron geometry. Models for the dip 
need extensions of the potential in Eq. ( PH ) , so that for the tetrahedron there are two energy 
scales present. Here several ways of achieving this goal are suggested: 

1. The effect of an isoscalar two-quark potential. 

As discussed in Ref. an isoscalar potential tu^ can be introduced into Vy - 
ensuring Vij = % for a colour singlet two-quark system - by the form 

V ij = ^r i .r j {v ij -w ij ) +Wij . (28) 

In this case, di — d% — since all of the are now equal to w and results in 
Ei = E 2 = E 3 = Aw. Therefore, from Table [IV] it is seen that the w range from 
-0.0035 to -0.0070 i.e. they have values much smaller than the corresponding 
= v given in the last column. A similar feature was found in Ref. jf], when 
the form in Eq. (ESI) was introduced to improve the model fit for squares and 



rectangles. However, as shown in Ref. [17|, in perturbation theory all terms of 

23 



0(a 2 ) are included in the two state model of Eq. (|T9|) with / = 1. Therefore, in 
the weak coupling limit Wij must be of 0(a 3 ) at least. 

2. The effect of a three- or four-body potential. 

The / factor is itself a four-body operator. However, it is conceivable that ad- 
ditional multiquark effects arise. Some perturbative possibilities are discussed in 
Ref. Jl7j. There it is shown that all three-quark terms arising from three gluon 
vertices always vanish, but that the four gluon vertex can contribute to 2-, 3- and 
4-quark terms at 0(a 3 ). However, in the tetrahedral case (r = d), cancellations 
result in this particular 4-quark term also vanishing. 

3. The effect of non-interacting three gluon exchange processes. 

These are also discussed qualitatively in Ref. [17] and contribute at 0(a 3 ) to 2-, 3- 
and 4-quark potentials. 

4. The effect of two quark potentials in which the gluon field is excited. 

The first excited state [^i*(r)] of the two-quark potential V\(r) is approximately 
given by V{(r) »s Vi(r) +ir/r - see for example refs. [§,0. Therefore, if a fourth 
state - based on such an excited state - is introduced into the model, since it is 
higher in energy than the three degenerate basis states so far considered, it will 
give attraction in the ground state. Furthermore, as the size of the tetrahedron 
increases this fourth state will approach the other three states, so that the at- 
traction felt in the ground state will increase - a trend seen in the tetrahedron 
results for E^ 2 in Table [n| 

The above isoscalar potential option now offers a reason for E\ = E 2 ^ 0. But, unfor- 
tunately, E 3 is still equal to E^ 2 since d\ = d 2 . However, there is no reason to expect any 
three or four body forces to also be purely isoscalars. In this case, their contributions to 
di and d 2 could be different and through the presence of the 1 — / factor in Eq. (|25|) any 
estimates of E 3 could be very model dependent. 
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VI. CONCLUSIONS 



This paper discusses three separate aspects of four-quark energies: 

1) Four-quark configurations (tetrahedra) involving three degenerate basis states. 
This showed the interesting result that - for tetrahedra - the ground and first excited states 
are degenerate in energy. The onset of this degeneracy can be seen in Fig. EL where - as a 
square (r = 0) gets deformed into a tetrahedron (r = d — 3) - the ground and first excited 
states, originally at -0.053 and 0.116, become degenerate at -0.026. This is a new feature 
not observed in earlier four-quark configurations (e.g. squares, linear,. . .). 

2) The use of the correlation matrix in euclidean time for extracting energies from the 
basic lattice data. 

Because the four-quark binding energy is the difference of two quantities that are comparable 
in magnitude, the final result involves a rather delicate cancellation - see Eq. flU). It was, 
therefore, considered worth while to investigate the effect of correlations in the basic data. 
Here it is shown that such an improved analysis of the lattice data leads to results that are 
essentially the same as those in earlier, less complete, analyses. 

3) The construction of a model in an attempt to explain these energies. 

Here it is seen that the /-model, introduced in earlier papers, needs to be modified before it 
can be used to described the tetrahedron results. Hopefully, a guide to these modifications 
can be suggested by perturbation theory, which has already proven useful in justifying the 
basic model at small interquark distances. 
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FIGURES 





FIG. 1. The four-quark geometry based on a square of side d parallel to the yz-plane and a 
distance r from that plane. A, B and C are the three possible partitions. The case when r = d is 
called the tetrahedron. 



FIG. 2. The binding energies - in units of the lattice spacing - of the four-quark states for the 
geometry of Fig. [l] for d = 3 and r = 0, 1, 2, 3, 4. 
Solid lines show lattice results: 
o - the ground-state binding energy E\. 
x - the first excited state energy E2. 

Dashed and dotted lines show model results with / = 1 from Eq. (|27j) for E\ and E2 respectively. 
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TABLES 

TABLE I. The binding energies Ek(d = 3a, r = 4a) obtained with method (I) from Eqs. 
ku is the number of energies extracted. 
T is the time range analysed. 
tim is the number of eigenvalues initially retained. 

nt is the actual number of unchanged eigenvalues finally retained using Um- 

X 2 /DoF is the ratio of Chi-squared/Degrees of Freedom with the actual result shown in the paren- 
theses. 

Case 8 is a fit with a unit correlation matrix and is equivalent to the usual Least Squares fit. 
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TABLE II. The binding energies Ei of four quarks in the geometry of Fig. || for d = 1 with 
r = 1, 2 and d = 2, . . . , 5 with r = d, d ± 1. Here r = d are the tetrahedra. 

1) ^4 + -B + C (I) is from Eqs. ©-(1) with A, 5, C shown in Fig. |. 

2) A + 5 + C (II) is from Eq. @ in a 3*3 basis - see Ref. §. 

3) A + B and -B + C are from Eq. (^) using only 2*2 bases 

4) V\ are the two-quark potentials i) V\{d,d) and ii) Vi(d, r) 

5) The symbol S indicates that the entry in the table is the same as that to its left - within 
numerical accuracy i.e. rounding errors. 

These results are from 2208 measurements contained in 69 blocks of 32 measurements each. 
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TABLE III. The notation is the same as Table || but with r = 1 and d = 2, 3, 4, 5 
These results are from 3008 measurements contained in 47 blocks of 64 measurements each. 
The Ei{d x d) are the corresponding energies for the nearby square of side d. 



(d/a,r/a) 


A+B+C (I) 


A+B+C (II) 


Ei(d x d) 


Vi 


(2,1) 




-0.0447(6) 


-0.043(2) 


-0.0588(3) 


i) 0.6689(4) 




E-2 


0.085(1) 


0.085(1) 


0.1414(8) 


ii) 0.6021(3) 




E 3 


0.69(3) 


0.66(3) 






(3,1) 


Ei 


-0.052(2) 


-0.049(2) 


-0.0531(5) 


i) 0.7974(8) 




E2 


0.089(3) 


0.09(1) 


0.1157(12) 


ii) 0.6992(5) 




£3 


0.53(2) 


0.55(5) 






(4,1) 


Ei 


-0.050(4) 


-0.047(1) 


-0.0524(10) 


i) 0.9102(15) 




E12 


0.088(6) 


0.088(3) 


0.097(2) 


ii) 0.7841(10) 




£3 


0.54(5) 


0.54(3) 






(5,1) 


Ei 


-0.052(5) 


-0.044(5) 


-0.047(3) 


i) 1.017(2) 




E2 


0.073(8) 


0.073(1) 


0.075(3) 


ii) 0.8652(8) 




E 3 


0.51(9) 


0.55(4) 
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TABLE IV. The energies for a selection of squares and rectangles in the ranges d = 1, . . 
and r = 1, . . . , 5. 

The Ei are from a 2*2 basis using 800(1600) for the rectangles (squares). In contrast to Ref. 
these results utilize Eqs. (||)-(||) and not Eq. @. 
The other columns use only 64 measurements. 

A + B + C denotes the 3*3 basis and A + B, B + C the corresponding 2*2 bases. 

The energies are those at T = 3 and not from any plateau in T. 

The notation is that: 

A has links along the sides of length r. 

B has links along the sides of length d. 

C is the state defined by the diagonals. 

S{L) is the state(s) with the lowest unperturbed energy. 
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